Two-gap superconductivity in MgB2: clean or dirty? 
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A large number of experimental facts and theoretical arguments favor a two-gap model for super- 
conductivity in MgB2. However, this model predicts strong suppression of the critical temperature 
by interband impurity scattering and, presumably, a strong correlation between the critical temper- 
ature and the residual resistivity. No such correlation has been observed. We argue that this fact 
can be understood if the band disparity of the electronic structure is taken into account, not only 
in the superconducting state, but also in normal transport. 
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Most researchers ascribe the superconductivity in 
MgB2 [0 to the electron-phonon mechanism, enhanced 
by interband anisotropy of the order parameter ^. 
Interband anisotropy, as expressed by the two-gap model 
Q , offers a simple explanation of many anomalous ex- 
perimental findings, most importantly of tunneling and 
thermodynamic measurements But there is a strong 
argument against it: As illustrated in Fig. |l|, existing 
bulk samples of MgB2 have essentially the same criti- 
cal temperature although their residual resistivities, po, 
vary greatly, between 0.4 and 40 ^fl cm. Such a be- 
havior is expected for s-wave pairing (Anderson's theo- 
rem), but not when two gaps are present. In that case 
one expects Tc to fall with increasing po- Indeed, im- 
purity interband scattering (magnetic and nonmagnetic) 
with rate 7 inter suppresses two-band superconductivity 
as: ATc oc 7inter /{t^Tc) and it is tempting to assume 
that 7intra ^ 7intcr OC po • For a sample with po ~ 40 
pf2 cm it seems unlikely that 7intor can be smaller than 
ttTc. In fact, the body of experimental evidence (Fig. 
|l|) can be reconciled with the two-gap model only if 
7intcr ^ 7intia- Until this paradox is resolved, the two- 
gap model for superconductivity in MgB2 cannot be ac- 
cepted, despite much compelling evidence. Two further 
problems are: (a) The high-temperature slope of the re- 
sistivity is clearly correlated with the residual resistivity 
(violation of Matthiessen's rule) [^, and (b) the plasma 
frequency estimated from the measured infrared reflec- 
tivity is 5 times smaller than the calculated one |@, |[ ^■ 

In this letter we shall show that the paradox can be 
resolved to support the two-gap model. It turns out that 
due to the particular electronic structure of MgB2, the 
impurity scattering between the a- and 7r-bands is excep- 
tionally small. Thus, the large variation of the residual 
resistivities reflects primarily a large variation of the scat- 
tering rate inside the a- and the 7r-bands, while the in- 
terband CTTT-scattering plays no role in normal transport. 
In the superconducting state, the two different gaps in 
the a- and the 7r-bands are preserved even in dirty sam- 
ples due to the extreme weakness of the crTr-interband 
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FIG. 1: Critical temperature for samples of varying qual- 
ity as a function of the residual resistivity. The theo- 
retical curves are computed in the two-band model, ac- 
cording to Ref. j^, with different ratios: Fo-o- /A'^o- (0) : 
F^^ /N^ (0) : /N^ (0) . Filled symbols refer to 'high- 
quality jingles': dense wires (A)[[1C| and single crystals 



(♦,■) JllJ 

samples 




. Half-filled symbols refer to 'high- Tc, high-p' 
Open symbols refer to samples of intermediate 



quality (0,o, A,V,n) p 



impurity scattering. 

MgB2 has two tt and three u-bands (Fig. ||) formed by, 
respectively, the two B Pz and the three bond-orbitals per 
cell, or, more correctly, by the corresponding Wannier- 
like functions. A bond orbital is the bonding linear com- 
bination of the two B sp^-hybrids which are directed 
along a B-B bond. The attractive potential from the 
Mg^+ ions in the hollows between the hexagonal boron 
layers is felt much stronger by a -electron than by a 
bond-electron and, as a result, the 7r-band is pulled so 
far down in energy that ~0.17 holes are left at the top 
of the CT-band. The strong coupling of these holes to 
the optical bond-stretching modes is what drives 
the superconductivity. Since the top of the a-band is 
at k|| = {kx,ky) = and is doubly degenerate, the holes 
are distributed in an upper heavy and a lower light band. 

The basic reason why crTr-impurity scattering is small 
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FIG. 2: LMTO bandstructure of MgB2 along the TA-line 
and in the plane {kz = j^) between the FMK and ALH-planes, 
where the a and n bands (fat) hybridize most. The FM/AL- 
direction is along, and the FK/AH-direction is perpendicular 
to a B-B bond. The orbital characters of the heavy and light 
(T-bands are explained in the text. 6x2x1 supercell bands 
for Mgi2B24 and MgiiB24 are shown along the main folding- 
direction, FM. For MgiiB24, two extra electrons and protrone 
dded, and the nuclear charge of each Mg increased by 2/11, 
to preserve the band filling and electroneutrality. 



is that the a and 7r-bands are formed from different lo- 
cal orbitals, and therefore are orthogonal on the atomic 
scale, rather than merely on an intermediate scale be- 
cause of Bloch factors. Moreover, the layered struc- 
ture and the compactness of the B 2s and 2p orbitals 
makes the cTTr-disparity in MgB2 much stronger than, 
say, the sd-disparity in a transition metal, where the sd- 
hybridization gap is almost as large as the d-bandwidth. 

Specifically, since a Pz-orbital has odd-parity, and a 
bond-orbital has even parity with respect to the B-layer, 
the only route for crTr-hybridization is via interlayer 
hopping, from a px-orbital in one layer to a bond 
orbital in another layer. The corresponding hopping 
integral, tbz, is, essentially, the geometrical average of 
the integrals ^ 0.1 eV and tjr^ ^ 1 eV, responsible 
for the fcz-dispersions of the a and 7r-bands |Q, and 
therefore small. Two further factors limit crTr-coupling: 
The second is that, in its interaction with the nearest 
bond-orbitals in the next layer, the B p^-orbital picks 
up merely the axial projection, which is essentially 
the s-character, on the boron above (or below) it. 
Near the top of the cr-band, the linear combinations 
of the three bond orbitals are, however, such that the 
contributions from the B s-orbitals cancel, so that the 
top of the cr-band is purely B px,Py-\ik.e. Hence, the 
only source of B s-character is tails of B p-orbitals 
centered at other sites. It turns out that the wave- 
functions for the heavy and light holes (z/ — h, I) are: 
|a„k) cx ^Tb-(r + ^-T)-p.(r-r-T)]e*-T, 
where T are the lattice translations, ±t are the posi- 
tions of the two borons in the cell (i.e., in a bond), and 
Ph/i (r) is a B p-orbital directed transverse/longitudinal 
to the k||-vector. From this representation, illus- 
trated in Fig. H, it may be realized that the B 
s-character often vanishes completely, and that it 
generally vanishes proportional to kf^ for the heavy- 
holes, and proportional to k„ for the light holes. The 



third limiting factor is the matching of the phase, 
(y9, between the two p^-orbitals in a bond, |7rzp,k) cx 
Et [Pz (r + r - T) e^'^C'") TP. (r - r - T)] e*-T^ 
and the phase between the corresponding B s- 
characters arising from the antibonding combina- 
tion, pn (r + r) — pn (r — t) . In the nearest-neighbor 
orthogonal tight-binding model for the 7r-bands, 
ip (k||) = arg {l -|- e^^ '^ + e'*' ''*^*'-' } , where a and b are 
the primitive translations of the layer. 

Due to their even/odd parity, the a and 7r-bands can 
only hybridize when kz ^ integer. Even then, as seen 
in Fig. H, the 7r+-band neither hybridizes with the heavy 
cr-band when k,, is along a bond, nor with the light cr- 
band when k|| is perpendicular to a bond. As may be 
realized from the pictures of the cr-orbitals (Fig. ||), 
the crossing with the heavy band occurs because the 
B s-character of that band vanishes exactly along this 
k-line, and the crossing with the light band occurs be- 
cause, along that k-line, the B s-character is purely an- 
tibonding between two borons, whereas the 7r_(_-band is 
purely bonding (ip ~ 0) . The two crTr-gaps seen in the 
figure are 0.2-0.3 eV, i.e., the hybridization matrix ele- 
ments, I (crk|iJ|7rk) I, are merely a per cent of the cr and 
TT bandwidths! 

We now discuss impurity scattering and use [15| : 
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for the rate of scattering to band n' of an electron in 
band n, by a weak localized impurity potential, V (r) . 
Here, denotes the average over the Brillouin zone, 
e„k is the band energy with respect to the Fermi level, 
and A^(0) = En^"(0) = EnU^i^nk) is the density of 
states per spin and cell. Typical defects for MgB2 are Mg- 
vacancies and Mg-substitutional impurities, which form 
easily, and B-site substitutions like N and C, which have 
a higher energy cost. The potential V (r) for a localized 
Mg-defect has the full point-symmetry of the site and, 
like the Mg^+ potential in the crystal, is felt more by a 
Pz-orbital than by a bond orbital. Hence, the largest 
matrix elements are those involving p^-orbitals near the 
impurity, i.e., the largest perturbation is of the energies 
of the p.-orbitals on the B hexagons immediately above 
and below the impurity, and of the corresponding t^^- 
This means that F^.^. should be large. Screening per- 
turbs the energies of the bond orbitals surrounding the 
impurity, and also perturbs t^, but to a lesser extent. 
Hence, we expect that Ttttt > Fo-o- for Mg-defects, albeit 
not for B-site substitutions. What contributes to T^tt, are 
matrix elements involving a pz and a bond-orbital, and 
most importantly, those on either side of a Mg-defect. 
Since this matrix element is the perturbation of tzb, it 
is expected to be intermediate between those of t^z ^^'^ 
t-^ij, like for the crTr- hybridization. Moreover, since the im- 
purity potential is fairly constant around a neighboring 
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boron, apz-orbital still picks up merely the B s-character 
which vanishes as kf^ for the heavy and as k„ for the light 
holes. This makes | (CTk|T^|7rk') | minute because kph and 
kpi are very small. Also the mismatch of phases between 
the a and 7r-functions will tend to reduce | (tTk|]/|7rk') |. 
Finally, squaring this small matrix element and inserting 
it in ([^), leads to an exceedingly small Fo-tt- 

To gain quantitative understanding of the disparity be- 
tween the scattering rates we have performed LMTO su- 
percell calculations for various impurities. Since the in- 
duced CTTT-gaps, 2| (fTk|y|7rk') |, are sensitive to their posi- 
tion within the cr-band (the B s-factor), we must choose 
a supercell which provides band- foldings near ep. The 
results shown in Fig. || were obtained with a 6 x 2 x 1 
supercell. The bands labeled Mgi2B24 are the same as 
those in the left panel, but folded into the smaller zone. 
The heavy c-band now crosses itself closely below ep, 
while the heavy-light and light-light crossings are a bit 
further down. The 7r_-band (fat) slightly above the top 
of the cr-band was originally at ML/2 and has been folded 
3 times into F. The MgiiB24 bands illustrate the effects 
of a Mg-vacancy: While the three 7r-bands get split by 
0.35 eV, and the heavy and light cr-bands by 0.27 eV (but 
by 0.04 eV at F), the crTr-splitting of the heavy band is 
merely 0.015 eV and that of the light band is merely 0.030 
eV! The squares of these splittings give estimates for the 
corresponding F's. For Mg-vacancies therefore, 

r.. >F,, »F,,. (2) 

We found very similar results for systems in which the 
Mg-vacancy was compensated by substitution of B by 
two C or one N: For MgisBsiN, the TTTr-splitting was 0.4 
eV, the CTcr-splitting 0.3 eV, and the crTr-splittings less 
than 0.03 eV. 

Let us now investigate how the relation (^ infuences 
the transport properties. These depend both on the 
impurity scattering and on the electron-phonon interac- 
tion (EPI). The interband anisotropy should be taken 
into account both in the impurity scattering (as outlined 
above), and in the EPI. The latter can be characterized 
by two sets of four spectral functions each: the stan- 
dard Eliashberg functions Q;^F„„/(a;), which define the 
superconducting properties and thermodynamical prop- 
erties like the electronic specific heat and the de Haas- van 
Alphen mass renormalizations, and the transport Eliash- 
berg functions al^Fnn' {'-^) ■ Of the calculated q;^F„„/ (w) 
functions (the details of the calculations are as in 
Ref. (14)), Q;^Fo-cr(w) exhibits a large peak at w « 70 
meV. Defining A„„' = 2 J w^^a^i^„„'(w)dw, we obtain 
the partial EPI constants, shown in Table |, which are 
similar to those obtained in In the following we as- 
sume that Pinter = 0, so the clean limit is appropriate 
(Tc is independent of the intraband F's). The super- 
conducting properties in the clean two-band model have 
been investigated in detail |l6j. Therefore we shall 
now concentrate on the normal transport. 
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FIG. 3: The DC resistivity in the clean (a) and dirty (b) case 
compared to experimental data for dense wires and for c- 
oriented films |7|, respectively. The lines are calculated in the 
effective two-band model with the indicated scattering and ab 
initio plasma frequencies Lof'— 4.14 eV, Lof'= 5.89 eV, u% = 
0.68 eV, and uj%= 6.85 eV. 

The explicit expression for the conductivity in the two- 
band model is |l^ (omitting Cartesian indices) 

l/PDc(T) = ^ ^ c.2;„/M/„(0,r), (3) 

n— cj,7r 

X [al{uj)F^^{uj) + ai.{uj)F^^{uj)] , 

where ja- = 7<T<T + 7o-7r, Itt = liTTi + l-na, and 7„„' ~ 2F„„/. 
Eq. (||) is essentially the standard parallel-conductor for- 
mula. Our assumptions are that 7^^. « and that ^„ 
differs much less than 7^ between 'good' (e.g., Ref.p^) 
and 'bad' samples (e.g. |^). 

The role of interband anisotropy, clearly visible in Ta- 
ble |, is different in superconductivity and electric trans- 
port. For instance, the critical temperature is given by 
the maximum eigenvalue of the A-matrix |l^, [l^ (i.e., 
mostly by its maximum element), while the conductiv- 
ity (H) is the sum of the partial conductivities. At high 
temperature, therefore, the slope dp/dT is determined by 

Y^n^lln/^trn, *-e., by the smallest Atr„ = Y.n' ^tmn'- 

In Fig.^ we show the temperature dependence of the 
DC resistivity for (a) a clean case with — Jtt ~ meV, 
and (b) a dirty case with 7^ = 54 meV and 7^ = 1.2 eV. 
In the two cases, all plasma frequencies are the same. 
The model is seen to describe both cases well. Note that 
7ct and 7^ determine not only the residual resistivity, but 
also the temperature dependence of the resistivity. In a 
one-band model, it would be impossible to reconcile the 
data of Refs. @, | with those of Ref. if they differ 
only by impurity concentrations, and the corresponding 
violation of Matthiessen's rule would be totally inexpli- 
cable. 

Why is the temperature dependence of the resistiv- 
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transport 


0.80 


0.41 


0.30 


0.15 


superconducting 


1.02 


0.45 


0.21 


0.16 



TABLE I: Superconducting and transport coupling constants 
A for the effective two-band model. The partial densities of 
states at the Fermi level for the two bands have values of 
iVa(0)=0.15 and iV^(0)=0.21 (states/ceU- spin-eV). 



ity so different in these two cases? Let us compare the 
clean limit, = In — 0, with the dirty-Mg-layer's 
limit, = oo, 7cr = 0. Of the two parallel conducting 
channels, in the former case the 7r-bands are responsible 
for conductivity at high temperatures, as was mentioned 
above, and even at T r:! the conductivity is mostly due 
to the TT-bands, their plasma frequency being higher than 
that one of the cr-bands. Since the EPI constant for the 
TT-bands is small, the temperature dependence of the re- 
sistivity is weak. On the contrary, in the dirty case, the 
TT-bands do not conduct, due to an overwhelming impu- 
rity scattering, and the electric current is carried only by 
the cr-bands. It is the strong EPI for this band which 
causes the temperature dependence of the resistivity in 
dirty samples.. 

To conclude, we suggest a new model for electric trans- 
port in MgB2. The main ingredients of the model are (i) 
interband impurity scattering in MgB2 is small, even in 
low-quality samples; (ii) intraband impurity scattering 
in the cr-band is small relative to the intraband TT-band 
scattering; (iii) high-resistivity samples differ from good 
samples mostly by the intraband 7r-band scattering rate. 
Of course, (iv) the phonon scattering is stronger in the a- 
band. This model explains well such seemingly inexplica- 
ble experimental facts as (1) absence of direct correlation 
between the residual resistivity and the critical temper- 
ature, expected in the two-gap model and (2) a strong 
correlation between the residual resistivity and the slope 
dp/dT in the normal state. Finally, we would like to 
point out that the existence of two qualitatively different 
scattering rates in the two bands should manifest itself in 
oth er cxpcrimcnta, auch aa optical and microwave apcc — 
tros copy, or Hall effect. In particular, accmingly myateri 
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